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Abstract 

We study collective excitations of cold (2+l)-dimensional fundamental matter living on 
a defect of the four-dimensional Af = 4 super Yang-Mills theory in the Higgs branch. 
This system is realized holographically as a D3-D5 brane intersection, in which the 
D5-brane is treated as a probe with a non-zero gauge flux across the internal part of 
its worldvolume. We study the holographic zero sound mode in the collisionless regime 
at low temperature and find a simple analytic result for its dispersion relation. We 
also find the diffusion constant of the system in the hydrodynamic regime at higher 
temperature. In both cases we study the dependence on the flux parameter which 
determines the amount of Higgs symmetry breaking. We also discuss the anyonization 
of this construction. 
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1 Introduction 


The gauge/gravity holographic duality has been recently employed to study compressible 
states of cold matter [I]^ . In these studies the intersection of two different types of branes 
(Dp and Dg with g > p) is considered. The higher dimensional Dg-branes are treated as 
probes in the gravitational background generated by the lower dimensional Dp-branes. On 
the held theory side the probe branes add hypermultiplets in the fundamental representation 
of the gauge group . These matter helds live generically in a defect of the unhavored theory. 
The probe approximation corresponds to the quenched approximation on the held theory 
side (see for an unifying formalism of the diherent brane intersections and for a complete 
list of references). 


In the brane setup the non-zero charge density needed to have a compressible state is 
generated by turning on a worldvolume gauge held [^. The dominant collective excitation 
of these systems at sufficiently low temperatures is a sound mode (the holographic zero 
sound [^). At enough high temperature thermal ehects dominate over quantum ehects and 
the system enters in a hydrodynamic regime in which a dihusion mode dominates. These 
two regimes are connected by a collisionless/hydrodynamic crossover transition. 


In this paper we will consider the intersection of D3- and D5-branes, according to the 
array: 


0123456789 


D3 : X 
D5 : X 


xxx______ 

XX _xxx___ 


This D3-D5 system is dual to a defect theory in which A7 = 4, d = 4 super Yang-Mills 
theory in the bulk is coupled to A/" = 4, d = 3 fundamental hypermultiplets localized at the 
defect [9,10 . We will restrict ourselves to the conhguration in which the D3- and D5-branes 


are not separated in 789 directions, which corresponds to having massless hypermultiplet 
helds. 


Turning on a hux of the worldvolume gauge held along the internal directions 456, one 

in which some components of the fundamental 


realizes the Higgs branch of the theory 11 


hypermultiplets acquire a non-vanishing vacuum expectation value. The worldvolume hux 
induces a bending of the D5-brane along the 3 direction. As shown in 11 one can then regard 
the probe D5-brane as a bound state of D3-branes or, equivalently, one can interpret that 
some of the D3-branes end on a D5-brane and recombine with it. The (2-|-l)-dimensional 
defect induced by the D5-brane represents a domain wall separating two regions with gauge 
groups with diherent ranks (the jump in the rank as one crosses the wall is proportional to 
the worldvolume hux). 


In this paper we study the collective behavior of cold matter conhned to the defect, when 
the system is in the Higgs branch. The hrst step in our analysis will be determining the 
precise conhguration of the probe which represents the system at non-zero chemical potential, 
temperature T, and magnetic held B (the T = B = 0 case was studied in [^). We will 
then study the spectrum of excitations and we will determine the dispersion relation of the 
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zero sound mode at T = 0. We will find a simple analytical expression for the speed of zero 
sound as a function of the flux. Moreover, since the intersection is (2+l)-dimensional, we 
can consider mixed Dirichlet-Neumann boundary conditions in the UV, which corresponds to 
performing an alternative quantization 13 and thus the charge carriers become anyons. In 


the presence of the magnetic field the spectrum of the zero sound mode is generically gapped 
although, as in [^, one can adjust the anyon parameter to some critical value such that the 
resulting spectrum is gapless. We will also study the system at non-zero temperature and 
we will find the corresponding diffusion constant. 


The rest of this paper is organized as follows. In section]^ we will determine precisely our 
brane configuration. The fluctuations of the D5-brane probe will be analyzed in section]^ In 
section we obtain the spectrum of the zero sound. Section is devoted to the calculation 
of the diffusion constant. Finally, in section we summarize our results and discuss some 
extensions of our work. 


2 The brane setup 


Let us consider the supergravity solution corresponding to a stack of N D3-branes at non¬ 
zero temperature. The corresponding near-horizon geometry is a black hole in AdS^ x S®, 
whose metric is: 

dslo = ^ ( - ^ , ( 2 . 1 ) 

where x = (x, ?/, z), = AnggNa''^ is the AdS^ radius and the blackening factor 


= 1 - d 


( 2 , 2 ) 


In (2.2) Th is the horizon radius, related to the black hole temperature T as = ttT. The 


D3-brane background is endowed with a Ramond-Ramond five-form whose potential 
will be denoted by The component of along the Minkowski coordinates is given 

by: 


^(4) 


- t^X 


R4 


(2.3) 


Let us now embed a D5-brane probe in the geometry (2.1) in such a way that it is extended 
along {x,y,r) and wraps a maximal C (parameterized by two angles 6 and (p). If the 
D5-brane is bent along the third Cartesian coordinate 2 ; = z{r), the induced metric on the 
worldvolume of the D5-brane is: 


dsl = r‘^[—fdt^ + dx^ + dy^] + 


r 1 


/ 


-f z'"^ 


dr"^ + dO"^ -|- sin^ 9 dp^ 


(2.4) 


where we have taken units in which the AdS^ radius R = 1. We switch on a worldvolume 
gauge field given by: 


F = A[dr A dt -|- B dx A dy + q sin 9 d9 A dp , 


(2.5) 
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with q being a constant (the amount of flux)j^ As usual, the rt component of F in (2.5) is 
required in order to have a non-vanishing charge density. The action of a D5-brane probe 
in the background geometry is given by the sum of the Dirac-Born-Infeld (DBI) and Wess- 
Zumino (WZ) terms: 


5'd5 — 


ifi y- det (9 + F) + Ts <iy (?<“' A F , 


( 2 . 6 ) 


where T 5 is the tension of the D5-brane and g is the induced metric on the worldvolume. 
The DBI determinant for our ansatz is: 

\/1 + f z'‘^ — \/l + sin 6 ^ , (2.7) 

while the WZ Lagrangian density is given by: 

Bwz = T 5 (74 A F = T^qr^ z' sin 9 dt A dx A dy A dr A d9 A d(p . (2.8) 


Therefore, after integrating over the angular variables, the total Lagrangian density becomes: 


C = dvrTs 

The equation of motion for At leads to: 

A', 


\/r^ + B‘^ a/I -|- — A't \/l + q^ + qr‘^ z' 


+ r'^ f z''^ — A't 


yi + y ^ d, 


(2.9) 


( 2 . 10 ) 


where d is a constant proportional to the charge density. From this equation we get A( as a 
function of z'\ 

- d^l + r^fz'^ 

* ^Jd? + (1 + g2) (r 4 + 52) ■ 

The equation of motion for z leads to the equation: 

f z' 


Vr^ + 52 


= ^/l + q‘^ + , 


( 2 . 12 ) 


a/1 -I- / 2:'2 — At 

where is a constant of integration. By imposing regularity at the horizon of the embedding 


function z{r) [1^, yields = qrj^. It is then possible to use (2.11) and (2.12) to obtain A( 
and z' as functions of the coordinate r: 


A' = ^ 

* r‘^ + d? + q^ + {1 + q^) B'^ 


z = 


q 


+ d'^ + q‘^rl + {1 + q^) B'^ 


(2.13) 


In what follows it is convenient to express the different results in terms of the chemical 
potential y ai T = B = 0, which is given by: 


4r(-)^ 

fi = At{r = oo)= / dr dK (2.14) 

do V ^ 

^The flux number q must satisfy the following quantization condition, q = Tia'k, with fc S Z (see, for 
example, [^). However, in units in which R = 1, the Regge slope is a' = l^/4^^Ngs. Accordingly, we will 
consider q as a continuous parameter. 


4 

































3 Fluctuations 


Let us allow fluctuations of the gauge field along the Minkowski directions of the intersection, 
in the form: 

A = + a{r,x^) , (3.1) 

where is the gauge potential for the held strength (2.5) and a{r,x^) = ay{r^x^) dx'^ is a 
huctuation. The total gauge held strength is: 


p _ -pA) _i_ 

tab - + 


(3.2) 


with is the two-form written in (2.5). In order to write the Lagrangian for the 

huctuations at second order, let us split the inverse of the matrix as: 


(()(“) + f<“>) ‘ = + j, 


(3.3) 


where Q ^ is the symmetric part and JL is the antisymmetric part {Q is the so-called open 
string metric). Then, the Lagrangian density for the huctuations is: 

r^ + B^ /_.. 1 


C 


+ {1 + q^) 


( gbd _ jac Jbd ^ i jcd jab^ ^ ^3 4) 


where the Latin indices take values in a,b,c G {t,x,y,r}. Notice that we are choosing a 
gauge in which Ur = 0. The equation of motion for a'^ derived from (3.4) is: 


da 


r^ + B^ 


+ d"^ + q'^ rf^ + {1 + q"^) B'^ 


^ca ~\~ ~ 3^^ 3^ 


Tca rrdb 


red /T<ih 


fab 


= 0 . (3.5) 


Let us write these equations in the case in which the huctuation helds only depend on 
the coordinates r, t, and x. We hrst Fourier transform the gauge held to momentum space 
as: 


at,{r,t,x) = 


duj dk 

(27r)2 


ttyir, uj, k) e 


—iujt+ikx 


(3.6) 


In what follows it will be understood that the gauge held is written in momentum space. 
Moreover, we dehne the electric held E as the gauge-invariant combination: 


E = k at + 00 


(3.7) 


The equations of motion reduce to a set of two coupled equations for E and the transverse 
gauge held huctuation ay. The equation for the huctuation of the electric held E is given 
by: 


E" + dr log 


r4 f + d"^ + q"^ + {1 + q‘^) B'^ + (1 + 5'^) + -B^) j 


+ 52 


(1 -I- g2)(ci;2 _ f k‘^)r‘^ -I- [(1 -I- g 2 ) 52 _j_ ^2j ^2 


E' 


1 {I + q‘^){oo^ - f + [{I + qy B^ + d^] oo^ 


J.A p 

AiBd 


E 


-|- d2 -)- g2 ^4 J^2 

(1 -I- q^ioj"^ — f kyy + [(1 -I- g^) B'^ + d^] 
r {y + B^f ^^4 + d2 + g2 rl + y + q^B^ (^d^ + (1 + q^ {y + B^'^ 


(3.8) 
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The equation for the transverse fluctuation is: 

a" + dr log + + + + 0^) 

1 (1 + g 2 ) ( a ;2 _ J + [^1 + ^ 2 ) ^2 ^ ^ 2 ] ^2 


r^P 


AiBd 


r'^ + cP + + (1 + q^) 

E 


r (r4 + B^)f ^r^ + P + q^rl + {l + q^)B‘^ ' 


(3.9) 


Let us now see how one can eliminate the dependence of Vh on the equations of motion by 
performing appropriate rescalings. First of all, we dehne a new radial variable f = r/r/j. 
Then, one can check that Vh is eliminated from (3.8) and (3.9) by dehning the new rescaled 
(hatted) quantities as: 


oj = — 
Th 


I- ^ 

k = — , 

Th 


3 ± 

^ „2 ’ 




(3.10) 


The resulting equations are obtained from (3.8) and (3.9) by taking r/j = 1 and substituting 
all quantities by their hatted counterparts. Hatted variables are utilized in the numerical 
integration of eqs. (3.8) and (3.9). 


4 Zero sound 


Let us now study the system at zero temperature. First we study the 
(3.8) and (3.9) near the Poincare horizon r = 0. Assuming that B is 
equations of E and ay are given by the coupled system: 


equations of motion 
small {B P), the 


E" + 


a 


// 

y 


+ 


4^2 . a;2 ^ 

_ j^f I _ rp 

r[P + B"^) P 

4^2 , a;2 

_l_ ^ 


4i BP 
r{P + i?2) 

r{P + B'^) 


(4.1) 


This is the same system as in the D3-D5 case with zero flux of [^. We can readily write its 
solution in matrix form as: 



r 


iuj \ S 
r ) (jj 



(4.2) 


where ci and C 2 are integration constants and we have imposed infalling boundary conditions 
at the horizon. Next we take the limit of low frequency and momentum in such a way that 
k ~ 0{e). For small uj the solution (4.2) can be written as: 


ui 



(4.3) 
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We now take the low frequency limit first. One can verify that the equations decouple in 
this limit. Actually, the equation for E becomes: 


F" + <9 lo? + (1 + g') ^ n 

® (l + g2)(a;2_ + 0 - 

This equation can be readily integrated as: 

(1 + - k‘^)p^ + o;^ (P 


E{r) = - Ce 


dp 


[d? + (1 + q^) p^] 
where E'^^'> = E{r ^ oo). Actually, if we dehne the integrals: 


(4.4) 


(4.5) 


/Ci(r) = 
/C 2 (r) = 


P 


^Jp"^ + d? [d? + (1 + q^) p^] 
dp 

[d2 + (1 + g2) p4] ’ 


dp 


then, E{r) can be written as: 

E{r) = E^^'> - Ce {I + q‘^){u‘^ - ICi{r) + d^ K. 2 {r) 


(4.6) 


(4.7) 


Let us now expand this result near the horizon. We hrst expand the integrals /Ci(r) and 
/C 2 (r) near r = 0 as: 


/Ci(r) = + 0(r2) , 


/C2(r) = K:2-- + 0{r^) , 


(4.8) 


where /Ci = /Ci(r = 0) and IC 2 = A^ 2 (?" = 0). It is interesting to notice that the quantities 
Ki and K 2 are not independent. Indeed, they satisfy the relation: 


(l + g2)£, + d^lC2 = 


dp 


P 


'0 \/ p^ + d? d 


(4.9) 


Moreover, if we dehne the hux function E{q) as: 


J^id) = ~ = (1 + -q"^ 


then, near r = 0 we can write E{r) as: 


E = CE^^r + ^ , 

d d I 




(4.10) 


(4.11) 


It is worth stressing that the whole effect of the hux in (4.11) is equivalent to multiplying 
k"^ by the hux function E{q). 

The equation for for low frequency is: 


a” + dr log () a' = 0 


(4.12) 
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This equation can be integrated twice to give: 


ay{r) = - Cy 


dp _ (0) Cy / 1 1 5_ \ 

^ ^ r V2’4’4’ r^y’ 


+ ci 2)2 

with a'y'* = ay{r ^ oo). For small r the previous solution becomes: 


(4.13) 


“»('■)» 4°’ 




(4.14) 


Let us now match the two expressions we have found for E and ay in this double low frequency 
and near-horizon limit (eqs. (4.3), (4.11), and ( 4.14[ )). Looking at the terms linear in r we 
get hnd relations between the constants Ci, C 2 , c^, and Cy\ 


u 

^1 7 ■) 

d 


C2 ^7 Cy . 


.UJ 

d 


The identihcation of the constant terms yields the following matrix relation 




Uy 




iB 

d 



d ^ d 


(4.15) 


(4.16) 


Let us now require that our fluctuation modes satisfy the following mixed Dirichlet- 


Neumann boundary conditions at the UV 15,16 


lim 


nr 


fr, - iwrq = 0. 


(4.17) 


with n being some constant (the Dirichlet boundary conditions correspond to taking n = 0). 
As in the case with g = 0, these conditions are equivalent to: 


lim E = —in lim Tr ^ a' 1 , 

r- —r*—^ 


lim a„ = i 


n 


CU ^ — k ‘^ r—>-00 


lim \r‘^E'] . (4.18) 


The quantities on the left-hand-side of (4.18) are the UV values E^^'^ and In order to 


obtain the values of the right-hand-side of the two conditions in (4.18), notice that the radial 
derivatives of E and ay can be obtained from their low-frequency values (4.5) and (4.13): 


dE 

dr 


— k^) cev 


-2 


day 

dr 


,-2 


(4.19) 


From these expressions we can recast the boundary conditions (4.17) as relations between 
the constants E^^\ a^\ ce, and Cy. Indeed, let us dehne En^ ando^ as: 


^( 0 ) ^ ^( 0 ) ^,nCy , 




inCE . 


(4.20) 


Then, (4.18) is equivalent to the conditions: 


4 “’ = 47 = 0 


( 4 . 21 ) 






















Moreover, combining (4.16) and (4.20) we conclude that En'^ and a^)i are related to the 
constants ce and Cy by the following matrix equation: 


^^(o)n 






— in 

a 


-t^+in 


liO I /i 

d d 



(4.22) 


Furthermore, the non-trivial fulhllment of the condition (4.21) is equivalent to the vanishing 
of the determinant of the matrix in (4.22), which determines the dispersion relation satisfied 
by uj and k for the zero sound modes: 

2 

— 2i^u/ + i ^E{q)ujk^ — ^ E{q) + (5 — ndf = 0 . (4.23) 


Notice that the effect of the flux in (4.23) is encoded in the substitution —)■ E{q) k‘^. Let 
us now solve (4.23) for a; as a function of k for small values of (ca, k). At leading order oj is 
real and given by: 




’42 ' 


+ 


{B - ndf 

/i 2 


(4.24) 


It follows from the last term in (4.24) that the spectrum is generically gapped for non¬ 


vanishing B and n. However, it can be made gapless by adjusting the alternative quantization 
parameter n to the critical value Ucrit = B/d. This fact is illustrated in Fig. where we 
compare the numerical results to our analytic formula (4.24). Moreover, from the coefficient 


of the momentum in the right-hand side of (4.24) we can extract the dependence of the speed 


of zero sound Uq on the ffux q. Indeed, by analyzing the behavior of the ffux function (4.10), 


it is easy to conclude that mq = Tl/v^ when g = 0, whereas it approaches the maximal 


possible value uq = ±1 as q ^ oo. Moreover, solving (4.23) for u at next-to-leading order, 
we find the attenuation of the zero sound: 


Imca =- 


l-|-g^ 5 3 2\ 1 2 . {B — nd^ 




(4.25) 


In Fig. 1^ we compare the analytic results Re a; and Imo; with the numerics with varying flux 
q and find very good agreement. 


5 Diffusion constant 


Let us now consider the fluxed D3-D5 system at non-zero temperature. First, we analyze the 
equations of the fluctuations near the horizon r = r^- It is easy to verify that the equations 
decouple in this limit and that the equation for E near r = takes the form: 


E" 


+ (^— + Cl) E' + ( 

\r — Th / — Th) T — Th/ 


E = 0 , 


(5.1) 
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Re w 


Re w 




Figure 1: (left) The dispersions for the anyon fluid as the effective magnetic held is decreased 
to zero. The red points are the numerical values for g = 1 and n = Q,ncritl‘^)'iT-crit (top- 
down). (right) The dispersions for the anyon huid at criticality n = Ucrit with increasing hux 
g = 0,1, 2 (bottom-up). Both hgures are done at d = 10®, B = 3 ■ 10^. 




Figure 2: We display the dispersions of the lowest quasinormal mode at low temperature 
(zero sound) as functions of hux g. The analytic formulas (4.24) and (4.25) (blue curves) 
match extremely well with numerics (red points) for the real and imaginary parts of the 
frequency a), respectively, as long as the temperature is kept small enough; we use d = 10® 
and A: = 10 with B = n = 0. 


where r^, A, ci, and C 2 are constants. Equation (5.1) can then be solved in a Frobenius series 
around r = r^. After this near-horizon expansion we perform a low frequency expansion by 
considering k rsj 0{e), oj C>(e^). The coefficients A, ci, and C 2 at leading order in e are: 


A = 


UJ 


IQrl 


Cl = 4(1 -F g^) 


-|- (1 -|- q^){r\ + i?2) 


^ + • • • 


C 2 = 


1 + g^ 


rh 


4 d? + {l + q^){r\ + B^) 


+ ... . 


(5.2) 


It can be checked that near r 
approximated as: 


= r/i, at leading order in e, the electric held E{r) can be 
Enh[l + /3(r - Th) ] , (5.3) 
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where E^h is the value of E at the horizon and /? is a constant coefficient given by: 


^ ‘'"^2 + (1 + qi ) ( r ‘ + B 2 ) ■ 


(5.4) 


We now perform the limits in opposite order. For low frequency, the equation of motion for 
E can be written as: 


dr log 




E' = 0 . (5.5) 


. ^r^ + (P + q^ rl + {l + q^)B^ {d? + (1 + q^) (r^ + . 

This equation can be integrated as: 

E(r) = + CEX(r) , (5.6) 

where E'^^'> is the UV value of E, ce is an integration constant, and X(r) is the integral: 

/•oo pi ^ 


X(r) = / dp 


^ + d? + q^r\ + {1 + q^) B'^ [ + (l + qi) (p4 ^2) 


(5.7) 


We now expand E{r) in (5.6) near the horizon: 

K + B^) CE 


E{r) = E^^^+CEX{rt,) 


[d 2 + (1 + g 2 ) (^4 ^ ^ 2 )] 


{r -Th) + ... . (5.8) 


Let us now compare (5.3) and (5.8). From the comparison of the constant terms we arrive 
at: 

EW = E^h-CEX{rh) , (5.9) 

while matching the linear terms yields: 

[<P- + (1 + <)2) (rji + B2)]" 


Ce ^ -P 


l 2 2 

^ ,2\ 'h 


ri + B^ 


Enh — i (1+? ^ 4 I D2 
oj ri + B^ 


Thus, we can write E^^'> as: 

E^^^ = Enh 


[d^+(l+g^) (r^+i?^)] ^ Enh 
(5.10) 


1+4 444 0’++St (1+«=)] h(a) 


( 6 . 11 ) 


From the condition E^^'^ = 0 we get a dispersion relation of the type u = —iDk^, with the 
diffusion constant D: 

2 

D = {l + «=) ^/EukXbXX+X) Hn.) . ( 6 . 12 ) 


In terms of the hatted variables dehned in (3.10), the diffusive dispersion relation can be 
written as: 

X = —i Z) P , (5.13) 
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with t he rescaled diffusion constant D defined as D = r^D = ttTD. It is immediate from 
(5.12) to get the value of Z): 


'+ il + B^)il + q^) , . 

D = + - -J{d,B,q) , 

1 + B ^ 


(5.14) 


where J'{d, B, q) is the integral 

J{d,B,q) = / dx - - . (5.15) 

+ g2 + (1 + g2) 52 [J2 + (^4 + 52)(1 + ^2)] 

In Fig. I^we compare the numerical and analytical results for D as a function of q for different 
values of the magnetic field B. 


Logb 



Figure 3: The numerical results (red points) for the diffusion constant H as a function of 
internal flux q follow the analytic results (5.14) (plotted in blue curves) spot on. The different 
curves correspond to varying magnetic field strength B = 0,1, 2 (top-down) with d = 10®. 


Let us study some limits of the formulas for the diffusion constant we have just found. 
First of all, it is interesting to point out that the integral J can be performed analytically 
when d = 0. The resulting expression for D is: 




1 + B2 V2’4’4 


(5.16) 


As d —)■ 0 when T —)■ cx), the high temperature limit readily follows from (|5.16|): 

.115 
'2’ 4’ 4 


D = v' 1 + 9 " F{1, 1 1 -q‘) 


Thus, at large T the diffusion constant behaves as: 




ttT 


T —)■ CX) 


(5.17) 


(5.18) 
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Interestingly, this last expression coincides with the one fonnd in 
the opposite limit T —)■ 0. We hnd: 


17 . One can also stndy 


p (1 + g2) (vrT)2 f cP 

2 {-kTY + B^ yd^ + {l + qYB^ 


i®-5- 

’ 4’ 2’ 


qY + 


2B^ 

IP 


T ^ 0 . 
(5.19) 


6 Conclusions and outlook 


In this paper we studied the collective excitations of cold holographic matter conhned to a 
(2+l)-dimensional defect of dd A7 = 4 super Yang-Mills theory in the Higgs branch. The 
string theory realization of the system is a D3-D5 intersection with flux on the worldvolume 
of the D5-brane. We found a simple analytic expression for the dispersion relation of the 
zero sound as a function of the flux (see eqs. (4.24) and ( 4.25[ )). The speed of zero sound 
and the attenuation grow monotonically as the flux increases. We also studied the diffusion 
constant at higher temperatures. 


Our work can be naturally extended along several directions. We could study other 
observables of the fluxed D3-D5 system with general boundary conditions. Some of these 
observables are the AC and DC conductivities of the anyonic fluid, as well as its diffusion con¬ 
stant. Moreover, we could easily extend our results to the general Dp-D(p-|- 2) intersections 
with flux. 


A more ambitious project could be the study collective excitations of the Higgs symmetry 
breaking even in more general holographic setups. One of such systems could be the D3-D7 
intersection with an instanton on the D7-brane worldvolume. The explicit expression of this 
instanton at zero temperature and non-zero density has been found in 12 . Moreover, it 
was argued in 18 that this setup realizes holographically the color-flavor locking phase of 


color superconductivity. The analysis of the collective excitations of this model is of obvious 
interest and we intend to address this problem in the near future. 
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